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Abstract. The possibility of the CDM axino is presented. From the early days of 
supersymmetry, axino has been considered as hot dark matter and warm dark matter. 
But the CDM axino has been considered quite recently. It is because the low reheating 
temperature of the universe is gaining momentum recently from various considerations 
in particle cosmology. I present a general introduction on the CDM's, in particular the 
axino, the calculation on the axino density in the universe, and the constraints imposed 
on the axino parameters. 



1 Introduction 

In this talk, I review the work collaborated with Covi et al. (!],§] on the possibility 
of the CDM axino. 

Modern cosmology needs dark matter and dark energy in the universe: ficnu — 
0.3, Qa — 0.65. There are several particle candidates for the CDM: heavy neu- 
trino in the GeV range, axion in the 10 -5 eV range, the lightest superpart- 
ncr(LSP) in the 100 GeV range, wimpzilla in the 10 12 GeV range, axino, and 
other hypothetical heavy particles. 

The particle mass and the interaction strength are provided by particle 
physics. In the standard Big Bang cosmology, this information leads to the de- 
coupling temperature (Td) below which the number of particles in the comoving 
volume is preserved. But the needed inflationary period reheats the universe, 
after a period of supercoolong, to the reheating temperature Tr. If Tjj > Tr, 
the particles are mostly diluted out and the first guess on its number density is 
that it cannot serve as dark matter(DM). But even after the reheating phase, 
the thermal production can be effective enough to create sufficient number of 
heavy particles to close the universe. It has been most extensively studied for the 
case of 0(100 GeV) gravitinoQ. Because of the gravitino problem, we consider 
the reheating temperature below 10 9 GeV.[] 

For Td <Tr, the DM candidates are 

• Neutrinos: Gf is the interaction strength. 

• LSP: -f-J- — is the interaction strength. 

For Tq > Tr, the DM candidates are 



1 For the wimpzilla case, the exit phase of the inflation produce the required density. 
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• Axion: y- is the interaction srength. Collective motion of axions below the 
QCD phase transion contributes. 

• Axino: -p- is the interaction strength. Thermal and nonthermal productions 
contribute. 



2 Axion and gravitino constraints 

Before discussing the axino CDM possibility, let us briefly review the cold axions, 
collectively moving in the universe. The basic difference of the axion candidate 
from the other ones is that axion is boson and it is possible to have a collective 
motion. The classical axion potential is extremely flat due to the very tiny invisi- 
ble axion mass[Q. Because of the extremely weak interaction strength, the axion 
vacuum stays at a point (a) for a long time. The vacuum starts to oscillate when 
the Hubble time (~ 1/H) is larger than the oscillation period (l/m a ),H < m a 
which occurs when the temperature is about 1 GeV(|]. We understand that the 
axion is the pseudo-scalar field, having the interaction only through the anomaly 

where F^ is the dual field strength (l/2)e tll/pa F f " J . We stress that there should 
not exist any other term in the axion potential. This kind of nonrenormalizable 
interaction can arise in several different ways: 

• From fundamental interaction such as in superstringQ: F a ~ Planck mass. 

• From composite models 0: F a ~ compositeness scale. 

• From renormalizable theories ||] : A global symmetry can have a gluon anomaly. 
If this global symmetry is spontaneously broken, there arises a Goldstone 
boson coupling to the gluon anomaly. F a ~ the global symmetry breaking 
scale. 

Then the current axion energy density is estimated to be|| 

Pa (T 7 ) = m a (T 7 )n a (T 7 ) ~ (^r^f • (2) 

If F a is large(> 10 12 GeV), the axion energy dominates the current energy 
density in the universe. Since the energy density is proportional to the number 
density, it behaves like a CDM. On the other hand, if F a is small then the axion 
interaction is relatively strong and too much axions are produced in the core of 
stars. SN1987A restricts its lower bound around 10 9 GeV, leading to the axion 
window 

10 9 GeV <F a < 10 12 GeV. (3) 



The gravitino produced thermally after inflation decays very late in the 
cosmic time scale(> I0 3 s), and can dissociate the light nuclei by its decay 
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products p0[. Inflation was used to dilute the primordial gravitinos, but thermal 
production of gravitinos are troublesome in cosmology if it exceeds || 

T R > 10 9 GeV. (4) 

Therefore, in supersymmetric theories we must consider relatively small reheat- 
ing temperature after inflation. 



3 Axino 

Therefore, in supersymmetric extension of the axion model, axion can be a CDM 
candidate but it is produced very late. Its supersymmetric partner is axino, and 
the reheating temperature Tr must be smaller than Tp. The axion a and axino 
a are gathered in a supermultiplet, 

<5 = —= is + ia) + V2a9 + F^OO (5) 
v2 

where s is the scalar partner of axion, the saxion. The axion couples to the gluon 
anomaly through (l/F a ){FF}. It is known that the potential arising from this 
gluon interaction settles (a) at zero, which is the Peccei-Quinn mechanism. 

For axino to be CDM, it must be stable or practically stable. Without the 
i?-parity conservation, this cannot happen. Thus, we require a practical R-parity 
conservation. For axino to be the LSP, it must be lighter than the lightest neu- 
tralino whose mass is expected to be around 100 GeV. Thus, an estimate of 
the axino mass is of prime importance for a CDM axino. Note that there is no 
theoretical upper bound on the axino mass. 

Since axion is almost massless, one might expect that the axino and saxion 
are almost massless in the first approximation. However, the saxion obtains a 
soft mass term below SUSY breaking scale. It is like the SM SUSY scalars. So 
the axino mass is intimately reated to the SUSY breaking scenario also. Actually, 
it is known that the axino obtains a substantial massjll]]. For a specific model, 
including the soft SUSY breaking terms, consider 

V = |/| 2 (|Si| 2 + \S 2 \ 2 )\Z\ 2 + (A x fS x S 2 Z - A 2 fF 2 a Z + h.c.) (6) 

where Z, Si, and 5*2 are the SM singlets, and the supersymmetric term comes 
from the superpotential W = fZ{S\S2 — F%). Since (Si) is of order F a , Z is 
of order the A term. Thus, the fermionic partners have the mass matrix of the 
form, 

51 / m a f F a \ 

5 2 ma fF a (7) 

Z V f F « f F a / 

where m a = f(Z). One eigenvalue is the axino mass ~ m„, and the other 
masses are of order F a . As seen from this example, the axino mass is basically 
a free parameter, but is expected to be somewhat smaller than the naive SUSY 
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breaking scale due to the coupling /p2[. But in some models, the axino mass 
can be much smaller than the SUSY breaking scale. Take a superpotential W = 
fZ(SiS 2 - X 2 ) + (A/3)(X - M) 3 with one more singlet X which carries a 
vanishing PQ charge. This superpotential is much more complicated to analyze, 
but still we can show that mg = 0(A— 2B + C) + 0(m 2 , 2 /F a ). For the standard 
pattern of soft terms, B = A — m 3 / 2 , C = A — 2m 3 / 2 - Then, the axino mass is of 
order keV. Thus, even the tree level axino mass needs the knowledge on the full 
superpotential ||l3f . Therefore, we can consider the axino mass a free parameter, 
and we restrict our discussion in the region where axino is the LSP, which is the 
most probable choice with the PQ symmetry. 

The universal axion coupling is to the anomaly term. Therefore, the ax- 
ino coupling to the gauge multiplet is the most important couplings. For the 
fcrmions, only the coupling to the top quark is important. We consider all these 
couplings when they become appropriate. 



4 Axino density in the universe 

The axino decoupling temperature is 

where Njjw is the domain wall number in axion models and F a n is the axion 
decay constant in units of 10 11 GeV. For Tn < Tr, at To the number den- 
sity is determined. For axinos not to close the universe, it should not exceed 
12.8 eV g*{To)/g e ff where = 915/4, g e ff — 1.5, which was used for the 
0(< 2 keV) warm DM axino by Rajagopal et aZ.|Q. But in our Tr < To case, 
due to the gravitino problem^, we do not consider this region seriously. On the 
other hand, we can consider much heavier axinos. The estimation of the axino 
density in the universe is considered for the thermal case and for the non-thermal 
case (by the neutralino decay). 

For the thermal production, we solve the Boltzmann equation, 

= y^ y W(i + 3 — ► a -I )v Te i)nirij + ^2(T{i — > a H ))rii. (9) 

The axino yield, the number density divided by the entropy s = {2 r K 2 / 45i)g s ^T , 
is split into two pieces one from the scattering process and the other by the 
decay, Y? p = ^ Yfj att + £\ Y % dec where 



yscatt = [ TR d (<T(»+j-^a + •••))"«", ydec = [** ,y (JV^ g + • ■ ■))? 

Jo sHT ' 1 J sHT 



For the scattering, we consider ten important processes shown in Table 1. On 
the other hand, the non-thermal production is basically by the decay process of 
the very long lived neutralino Q. These contributions are shown in Fig. 1. For a 
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Table 1. The cross sections for the different axino thermal production channels in- 
volving strong interactions. Masses are neglected except for the plasmon mass m e //. 
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low reheating temperature, the thermal production is not effective. In Fig. 1, the 
bulge at the lower left corner is the non-thermal production. In Fig. 1, we show 
for F a — 10 11 GeV so that the cold axions are not the dominant component of 
CDM. 

For the calculation of the non-thermal production, we consider LOSP and 
NLSP where LOSP is the lightest superpartner among ordinary SM particles and 
NLSP is the next-to-LSP or the second lightest superpartner. It is most likely 
that the LOSP is the lightest bino-like neutralino and NLSP can be LOSP or 
gravitino or something else. For a bino-like LOSP x, the decoupling temperature 
of the LOSP is T D ~ m x /20. At T < T D and r x < H, the LOSP decays during 
the radiation dominated era to produce axinos through \ — > 5, + 7 to yield Q 



Qa = ^-n x . (11) 



Therefore, we have a simple expression, 

m 

The neutralino decays at the cosmic time 
f n«q 1 A/128 \ 2 f FjNnw V A00GeV \ 



Depending on the models, the decay modes can be different. If there are more 
channels for the x decay, then it will give an even more efficient implementation 
of the non-themal production of axino through the neutralino decay. 



6 Jihn E. Kim 



5 Constraints 

We must impose the conditions that (1) not too much axino energy density at the 
time of Big Bang nucleosynthesis (BBN), (2) the BBN not spoiled by x decay 
producing the SM particles, and (3) axinos becoming cold before the galaxy 
formation era. 

Note in passing that if < 77i 3 / 2 < m x , i.e. for the gravitino NLSP then 
the gravitino problem disappears since it decays to axion and axinojl5). On the 
other hand, if \ is the NLSP, the gravitino problem is prese nt |3| . In this case, 
the non-thermally produced axinos are estimated to give Eq. (|ll|). In Fig. 2, we 
include all these constraints and plot the restriction on the reheating temperature 
versus axino mass. This figure is not changed even if there is no gravitino problem 
because of the gravitino NLSP. The shaded region is where the axino can be a 
CDM possibility. This happens O(GeV) axino. 

Note that for a high reheating temperature the thermal production contribu- 
tion dominates. In the figure, there is the solid line denoting the critical energy 
density by axinos. Even if the reheating temperature is below the critical en- 
ergy density line, there still exists a axino CDM possibility by non-thermally 
produced axinos |Q. Note that non-thermally produced axinos can be CDM for 
relatively low reheating temperature(< 10 TeV) for which axino mass should be 

10 MeV < rria < m x : non — thermal axinos as CDM possibility (12) 

The shaded region corresponds to the MSSM models with the constraint fl x h 2 < 
10 4 but still allowable axino energy density. One can choose 30 % CDM in this 
region. However, if we restrict further that all SUSY mass parameters are below 
1 TeV, then we have fi x h 2 < 10 2 . For a sufficient axino energy density for the 
axino CDM, then we have 

m a > 1 GeV (13) 
but less than the neutralino mass. 



6 Detection possibility of CDM axino 

If i?-parity is exact, then the CDM axino cannot decay, and there is no way to 
prove the CDM axino possibility. If there is any chance to prove the CDM axino 
by observation, then the axino must decay [jl6). For it to constitute most of CDM 
in the universe, its lifetime must be of order the age of the universe. Since it is 
assumed to be the LSP, the i?-parity must be feebly broken. Then the CDM 
axino can decay. 

For the simplicity of discussion, let us consider the bilinear terms only 

W = fi a L a H 2 (14) 

where a(= 1,2,3) is the flavor index, fi a — O(eV) is the constraint from the 
neutrino mass bound. This bound applies to the heaviest neutrino, presumably 
v T . In this case /U3 is bounded as 

iMsl < M l J 2 MeV. (15) 

H2,TeV r 
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With the i?-parity violation, we expect the decays 

a — > v + 7(01' l + l~), 



a — > v + a (16) 

a — > T + + 7T~. 

In fact, the 1/7 mode is most important in cosmology. For the cosmic impor- 
tance of the decay, we need to know the lifetime. Let us try phenomenological 
Lagrangians, 

e o 0oV>, ie 1 -^F' J ' l/ aj 5 [j fl ,'y u ]ip ) (17) 

fa 

for the final pseudo-scalar and photon, respectively. Then, for each case the 
lifetime can be estimated. In particular, to the photon mode 

4xl0 4 o , lrl 

n = ^s] (18) 

where Pi is the phase space factor. e± is known if the i?-parity violation is given. 
Since the most stringent bound comes from the diffuse gamma ray, we focus on 
this mode. The observed gamma ray flux is 



< (10- 3 - lO-^Vcm^sr-Ms- 1 ]. (19) 



If tree < T a < to where to is the age of the universe, most axinos have decayed 
by now and the flux has a peak at Eq = (m a /2)(T a /to) 2 ^ 3 ■ Since the photons 
from axino decay should not exceed the observed flux, we obtain 

r a < lQ- w t f2 3 J 2 h 3 (20) 

which is in conflict with r a > t rec ~ 10 13 s. 

If T a > to, a similar consideration gives the axino lifetime larger than 4 x 
10 24 ~ 26 Q a h s. This translates to the 113 bound, 

to < 10 2 ~ 3 eV. (21) 

If [13 in ( ^l| ) is taken, the idea of neutrino mass generation by i?-parity violation is 
not enough. Thus, the unstable CDM axino is not viable with the neutrino mass 
generation by i?-parity violation. Still, the CDM axino decay can be detected if 
the parameters lie in the region given above. 



7 Conclusion 

In conclusion, we discussed 

• Supersymmetry is the solution of the gauge hierarchy problem, and the PQ 
symmetry with a very light axion with 10 9 GeV < F a < 10 12 GeV is the 
solution of the strong CP problem. This leads to the axino which can be the 
LSP. 
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• The gravitino problem must be considered in this case, with Tr < 10 9 GcV. 

• Due to the much more stronger axino production than the gravitino produc- 
tion, Tr must be much smaller than the one given by the gravitino consid- 
eration. It can be in a multi TeV region. 

• Thermal and nonthermal production of axinos are possible. In this case, 
O(GeV) axino CDM possibility exists. 

• Finally, the attempt to explain the neutrino mass via the i?-parity violation 
mechanism is inconsistent with the CDM axino. However, in this case there 
is a room for the detection of axino decay debris. 
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Fig. 1. Y- a as a function of T R for F a = 10 11 GeV and m 9 - = 1 TeV. 
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Fig. 2. Constraints on Tr and axino mass. 



